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The nonlinear response of an adsorbed layer on a periodic substrate to an external force is studied 
via a two dimensional uniaxial Frenkel-Kontorova model. The nonequlibrium properties of the 
model are simulated by Brownian molecular dynamics. Dynamical phase transitions between pinned 
solid, sliding commensurate and incommensurate solids and hysteresis effects are found that are 
qualitatively similar to the results for a Lennard- Jones model, thus demonstrating the universal 
nature of these features. 



I. INTRODUCTION 

In recent years, there have been an increasing effort 
at trying to understand the frictional force between iwo 
macroscopic sliding surfaces at the microscopic levelErH. 
It is now understood that in the case of boundary lu- 
brication, the tightly bound monolayer in the regions of 
closest contact between the two surfaces plays a crucial 
role in determining the macroscopic sliding friction. At 
the same time, new experimental techniques such as the 
quartz-crystal micrpbalance (QCM)E2l and atomic force 
microscope (AFM)H allow the direct probing of frictional 
force at the atomic level. Thus, the central theoretical 
issue in understanding boundary lubrication and atomic 
scale friction is to determine how an adsorbed layer re- 
sponds to an external driving force. In the linear regime, 
the response is characterized by the collective diffusion 
constant which is itself an important transport coefficient 
in determining surface dynamics such as growth of thin 
filmsl-j. The nonlinear regime is the one relevant for a 



tual conditions of macroscopic sliding motion. Persso 
has made a series of pioneering studies on this problem 
based on a model of an adlayer of particles interacting 
with Lennard-Jones potentials. The central result is that 
there are intriguing dynamic phase transitions between 
pinned solid, sliding solid and liquid phases. Strong 
hysteresis effects exist in the transitions between these 
phases, leading to the 'stick and slip' motion often ob- 
served in the sliding motion between two macroscopic 
flat surfaces under an external driving spring. In addi- 
tion, Persson has used a general hydrodynamic argument 
to predict that the ratio of the static and sliding friction 
threshold should take a universal value of approximately 
2 for all systems with strong coupling to the substrate. 
The qualitative features of the results from this 



isotropic Lennard-Jones system are very encouraging. 
Besides the relevance for the boundary lubrication prob- 
lem, the dynamic phase transitions that have been ob- 
served are of intrinsic interest and very different from 
the more familiar equilibrium transitions. However, ac- 
tual adsorbate layers usually have strong anisotropies and 
complicated many-body interactions among the adatoms. 
The question then arises as to how universal are the 
results for the Lennard-Jones system. Recently, Braun 
et au have studied the nonlinear mobility of an adlayer 
based on a generalized vector Frenkel-Kontorova model. 
They have also observed a strong hysteresis effect in the 
transition between the pinned and running state, while 
the nature of the running state and the ratio of the static 
to sliding friction threshold they obtained are quite differ- 
ent from the corresponding results of Persson. However, 
it is difficult to make a direct comparison between these 
studies. Persson focused mainly on an initial pinned state 
that is commensurate with the substrate, while Braun et 
al focus almost exclusively on an initial state that is in- 
commensurate with the substrate. In the latter case, the 
sliding motion is dominated by the response of the kinks 
or antikinks to the external force. 

To facilitate a more direct comparison with Persson's 
study and address the question of universality in the 
salient features, we present in this paper the results from 
a study based on a uniaxial Frenkel-Kontorova model 
that is different from the one employed by Braun et al. 
In particular, the particles in our model have harmonic 
spring-like interactions that prevents them from moving 
away from each other. This implies a fixed coordina- 
tion number for each adparticle and a very different na- 
ture of the liquid phase from the corresponding phase for 
the Lennard-Jones system. Moreover, the particle dis- 
placement is a scalar variable in contrast to the vector 
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displacements in the other models. Also, in this paper 
we focus mainly on an initial pinned state that is com- 
mensurate with the substrate. Surprisingly, we find that 
almost all the features of the dynamic phase transitions 
and hysteresis effects found in the Lennard- Jones system 
are present also for our generalized Frenkel-Kontorova 
model. The only notable difference between the models 
is in the value of the ratio of the static to sliding friction 
threshold. We present below the details of the model and 
the results. Preliminary results have been reported in an 
earlier conference proceedings. 

II. MODEL 

We consider a two-dimensional version of the Frenkel- 
Kontorova (FK) model, which-is the simplest extension 
of the one-dimensional modelij, obtained by coupling 
linear harmonic chains along the additional dimension. 
The particle displacement is always along the same direc- 
tion which makes it particularly convenient for numerical 
computations. This model is clearly more appropriate to 
describe uniaxial systems of adsorbed layerstfl however 
we believe that the main features of the dynamical be- 
havior which we consider in this work should be indepen- 
dent of this type of anisotropy. We have in fact checked 
that the dynamical transitions and hysteresis effects re- 
mains essentially the same if a more isotropic version of 
the modelij is used. 

The model is described by the Hamiltonian 

+ ^(^ +1 -^) 2 -C/ cos(^)] (1) 

where m is the mass of the particles; U the amplitude 
of the periodic potential with period a ; K x and K y are 
elastic constants in the x and y direction and b is the 
average distance between the particles in the absence of 
the external potential. The second term in the brackets 
represents the elastic interaction within a chain and the 
third term the interaction between adatoms on different 
chains. Note that the index i, j labels the particles within 
a chain with integer i [i — 1, N x ) and the chain with 
integer j (j = 1, N y ), so that the system consists of 
a total number of N x x N y particles. In the simulations 
we have considered N x = N y = N and K x = K y — 
K . Typical equilibrium particle configurations above and 
below the depinning temperature T c are illustrated in 
Fig. [j]for b/a = 2 

We use periodic boundary conditions in most of the 
calculations. This is obtained by requiring that the par- 
ticle (N x + 1, j ) is a periodic image of the particle 
(1, j) in the same chain and so xm^+i.j = + N s a, 
where N s is the number of local minima of the peri- 
odic substrate in the x-direction. In the other direc- 
tion, periodic boundary condition means that x^jv +i 
= Xi i. The average distance between particles, b, in 



each chain is determined through the boundary condi- 
tion and is given by b — N s a/N and therefore the over- 
layer coverage, or atomic concentration, can be defined 
as 9 = N/N s = a/b. With this set up, 9 and b cannot be 
changed independently. With free boundary conditions 
in the ^-directions, the ratio b/a can assume arbitrary 
nonrational values. The system now can remain com- 
mensurate for a given. 9 only up to a critical value S c of 
the misfit parameteiO i5 = (b — a)/a. This boundary con- 
dition is particularly useful to study incommensurability 
effects. 

An important feature of the model is that the parti- 
cles have a fixed number of neighbors. So, even at high 
temperatures, where a liquid phase with particles diffus- 
ing around independently is usually expected, there is 
some crystalline order although the correlations decay as 
a power lawlifl. We shall, nevertheless, refer to this phase 
as a "fluid ", keeping in mind its distinct nature. At 
low temperatures, the periodic potential is able to pin 
the overlayer and true long-range order prevails leading 
to a finite static friction. The low and high tempera- 
ture phases are separated by a commensurate solid to 
fluid transition at a critical temperature T c where the 
overlayer at high temperatures become disordered and is 
effectively depinned from the substrate. We shall also 
refer to this transition as "melting" although topologi- 
cal defects, as dislocations, are not allowed due to the 
fixed- neighbor constraint. In the absence of a periodic 
potential, dismxLsr plays an important role in pinning 
the overlayenl3il3 and glassy behavior is expected. 

The periodic potential in Eq.(l) arises from an adia- 
batic average over the substrate degrees of freedom. The 
nonadiabatic excitations of the substrate provide energy 
exchange between the substrate and the adlayer and lead 
to damping and fluctuation in the motion of the adatoms. 
In our model, this effect is described by the usual Marko- 
vian Langevin dynamics. The substrate acts as a heat 
bath at an equilibrium temperature T and the coupling 
of the substrate excitations and adlayer is characterized 
by a friction parameter r\. Note that in the absence of the 
periodic potential, the adatom interactions do not affect 
the collective response of the adlayer and the macroscopic 
sliding friction is equal to the microscopic coupling pa- 
rameter rj. . 

III. GENERAL PROPERTIES OF SLIDING 
FRICTION 

To characterize the response of the adlayer to the exter- 
nal driving force, we define the sliding friction coefficient 
fj through the relation 

v d = Fj mfj (2) 

where Vd is the average velocity of all the adparticles rep- 
resenting the drift velocity in the direction of the force. In 
general, the sliding friction fj has a complicated depen- 
dence on the external force F in the nonlinear regime. 
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First, we remark on a number of qualitative features of 
the drift velocity Vd and the sliding friction fj that-,can 
be readily understood without detailed calculationsta In 
the presence of the external force F, the potential is tilted 
like a washboard with local minima which become weaker 
as the strength of the force is increased. Beyond a critical 
force value F Q , the effective potential minima disappears. 
For the present model, -Fq-^p 2%U , the same result as for 
the one-dimensional caseli-i At T = and for a com- 
mensurate value of coverage 9 , Vd remains zero for any 
applied force less than F a . When F >> F 0) the layer 
slides as a whole with a sliding friction fj > r\. At low 
but finite temperatures, the drift velocity induced by a 
small external force is limited by the nucleation and mo- 
bility of defects in the commensurate structure. In one- 
dimension, these defects are kinks and antikinks (local 
compression or extension of the commensurate atomic 
configuration) and can always be excited at any finite 
temperature leading to a finite fj at any T ^ . In two 
dimensions, however, the corresponding excitations are 
linear defects (domain walls) which have a nonzero free- 
energy per unit length for temperatures below the melt- 
ing transition temperature T c . As thermal fluctuations, 
these defects can only be nucleated as closed loops with 
an activation energy proportional to its length and there- 
fore do not contribute significantly to the linear mobility 
of the system, lim^^o Vd/F, in the thermodynamic limit 
and one expects that for a pinned solid, 1/fj is essentially 
zero for T < T c in this regime. When the external force 
is increased from zero, activated nucleation and growth 
of commensurate domains induced by the external force 
can contribute to a finite mobility, leading to a nonlin- 
ear friction fj > r) with a complicated dependence on the 
external force F. Eventually, in the limit of large exter- 
nal force such that the effect of the periodic potential is 
negligible altogether, fj again reduces to the microscopic 
friction parameter rj. 

In the nonlinear regime, the behavior of drift velocity 
Vd for increasing external fqtce F is expected to depend 
strongly on the initial phasea at F = 0. If the adsorbate, 
when F — 0, is a fluid, Vd will be nonzero for arbitrar- 
ily small external force. This result is also true for the 
present model. In the fluid phase above T c , the over- 
layer is essentially depinned and any small force leads to 
nonzero Vd- Moreover, the relation Vd — f{F) should 
show no hysteresis. For an initial state at F = that is 
incommensurate, the critical value of the force for finite 
mobility can be much lower .as shown in previous works 
of PerssonB and Braun et al.Q. 



adsorbates. The equation of motion for the particle with 
coordinate Xij is given by the Langevin equation 

9V dU „ , 
mxij + mrjXij = —- hf + (3) 

where U is the periodic potential, V is the adsorbate- 
adsorbate harmonic interaction potential and is a ran- 
dom force that is related to the microscopic friction pa- 
rameter rj by the fluctuation dissipation theorem 

< fij(t)fi>j>(t') >= 2 n mk B T S itV S jtf S (t - (4) 

We have studied the sliding friction by simulating the 
above equations using methods of Brownian molecular 
dynamicsE-D which treat accurately the effect of the sys- 
tematic and stochastic forces in the problem. The cal- 
culations were performed at a commensurate coverage 
8 = 1/2 (b/a = 2). Typically, the time variable was 
discretized with time step St = 0.02 — 0.06 r where 
t = (ma 2 /Uo) 1 / 2 and 4— 10 x 10 5 time steps were used in 
each calculation allowing 10 5 time steps for equilibration. 
Calculations were performed as function of temperature 
and external force on systems containing N x N adsor- 
bate particles, where typically N — 10, and the elastic 
constant was set to K = 10 . 

In the calculations, the system was allowed to evolve 
into a steady state such that the time average of physical 
quantities, like the drift velocity Vd, approached a con- 
stant. The effective temperature T* of the overlayer dur- 
ing sliding was obtained by two different methods. Using 
the fluctuation of particle velocities from the average as 

kT* = m(< v 2 > - < v > 2 ) (5) 

and by equating the energy transferred to the adsorbed 
system by the external focce to the energy transferred to 
the substrate which givesa 

T* F 2 n 

■^ = 1 + — (6) 

where fj = F/mvd- The two methods agree within the 
estimated errors, for calculations expected to be equili- 
brated in a steady state, and provide additional support 
for the assumption that local thermalization of the over- 
layer can be described by an effective temperature T*. 

V. RESULTS AND DISCUSSION 



IV. MOLECULAR DYNAMICS SIMULATION 

To analyse the nonlinear sliding friction of the model 
presented in section II, we have performed a molecular 
dynamics simulation study of the Langevin equation re- 
sulting from the model. An external driving force F, 
representing the effect of the other surface in the bound- 
ary lubrication problem, is assumed to act on each of the 



In this section, we present results from our molecular 
dynamics simulations. We use units in which a = 1, 
m = 1 and Uq = 1 which sets the scales for length, mass 
and energy respectively. The corresponding time scale 
is t = (ma 2 /Uo) 1 / 2 , and the temperature scale is To = 
Uo/ks- It is also convenient to normalize the external 
force by the zero-temperature critical value F a . 

To obtain an estimate of the critical temperature T c , 
we calculated fj as a function of temperature for a small 
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value of the external force, F/F a = 0.035. At low tem- 
peratures, l/fj is essentially zero and the adlayer is in a 
pinned state. As the temperature is increased 1/fj re- 
mains zero untilEl T c ~ 2.3, beyond which it increases 
significantly and should approach 1 at high enough T. 
The value of T c correspond to the commensurate solid- 
fluid transition. Above T c , the overlayer is effectively 
depinned from the substrate. 

In Fig. U we show the behavior of Vd and the effective 
temperature of the layer T* , obtained at a temperature 
T = 3 which is above T c . In this case, the adlayer is 
in a fluid state at zero external force. The results for 
F increasing from zero and decreasing from its highest 
value are the same, showing a total absence of hysteresis 
effects. The effective friction coefficient fj is finite even 
at vanishingly small values of F. In the absence of the 
periodic potential, fj = rj and Vd should depend linearly 
on F as indicated by the dotted line in Fig. |[ In the 
present case, the net motion of the overlayer is limited 
by defect excitations induced by the coupling to the pe- 
riodic potential leading to an enhancement of fj over the 
microscopic friction parameter rj and an increase of T* 
over the equilibrium temperature T, as F increases. The 
sliding friction fj approaches r\ only for larger F . The 
calculation of the effective temperature T* assumes that 
the particle velocities have a Maxwellian probability dis- 
tribution of width ksT* . In fact, as shown in Fig. the 
velocity distribution can be well described by a Gaussian 
distribution. 

For T < T c however, where a commensurate solid 
phase prevails when F — 0, a critical force F a is nec- 
essary to depin and initiate sliding. It was found for the 
Lennard-Jones systema that in this case the nonlinear 
sliding friction exhibit hysteresis as a function of F, i.e., 
the relation between velocity and external force depends 
on whether F increases from zero or decreases from a high 
value. We find the same behavior in our model for suffi- 
ciently low temperatures. In Fig. ^, we show the results 
for Vd and the effective temperature of the adlayer T* as 
a function of F for r\ = 0.6 obtained at a temperature 
T = 0.5 which is much less than T c . This temperature 
is also considerably below the amplitude of the periodic 
potential U a . As will be discussed below, we find that 
T < Uq seems to be necessary in order that a sharp hys- 
teresis feature appears in the Vd x F characteristics. The 
calculation was performed with the adlayer initially in 
a pinned state, then the external force was increased to 
different values. An hysteresis loop is clearly seen in Fig. 
|] and share the same features .as found in the simulations 
of the Lennard-Jones systemsH. 

In Fig. ||, it is shown that the initial sliding phase, 
when the applied force is increased beyond the static 
threshold F a , has almost the same temperature as the 
substrate and the substrate potential provides no addi- 
tional resistance. This is a sliding solid phase. It corre- 
sponds to the floating or incommensurate solid phase in 
the equilibrium situation. It is important to note that 
this phase is dynamically generated by the external field 
since the initial temperature is below T c and a commen- 



surate solid phase would be the equilibrium state. With 
the choice of periodic boundary conditions in our model, 
the misfit 5 is zero and the configuration of the sliding 
incommensurate solid phase in Fig. ^ is indistinguish- 
able from a commensurate phase. However, using free 
boundary conditions and varying the misfit parameter 
we can confirm that this phase corresponds to a sliding 
solid phase incommensurate with the substrate. In Fig. 
U we show the structure factor S(q) for two values of the 
overlayer lattice spacing, b/a — 2 and b/a = 1.9 , and 
at two different limits of the external force, F/F„ = 0.1 
in the pinned phase (Fig. |^b) and F/F a = 0.7 in the 
sliding phase (Fig. ||a). In the sliding phase, S(q) is 
peaked at the wavevector 2n/b independent of the sub- 
strate while it remains locked to the commensurate value 
7r in the pinned phase. Thus the external force induces 
a pinned commensurate solid to sliding incommensurate 
solid transition at F — F a . As the force is decreased from 
its highest value through F a , there is a jump in T* accom- 
panied by a simultaneous decrease in Vd at a critical force 
F c . This corresponds to an overheating of overlayer to 
an effective temperature which is above T c , dynamically 
"melting" the sliding solid phase. This can be regarded 
as a sliding solid to fluid transition. Finally at F = F},, 
the fluid phase condenses into a pinned commensurate 
solid phase and the temperature of the adlayer and the 
substrate become the same again. The fact that there 
is a condensation from the "melted" phase is confirmed 
by noting that the temperature of the adlayer at this 
condensation is almost identical to the equilibrium crit- 
ical temperature T c . Perssonfl has observed very similar 
dynamical transitions for Lennard-Jones systems, show- 
ing that these features are universal in the non-linear re- 
sponse of overlayers and not just for specific models. In 
fact, other versions of the Frepkel-Kontorova model also 
show some of these features!! In addition the present 
results show that a conventional "liquid" phase is not 
required and a high temperature thermally disordered 
phase with nonzero linear mobility is sufficient to show 
some of these effects. 

We now consider the behavior of the nonlinear sliding 
friction at a higher temperature T — 1.5 but still less than 
T c . As can be seen from Fig. |6[ the hysteresis disappears 
but F a can still be defined. In addition, a smooth sliding 
fluid to solid transition takes place at F c > F a . At this 
higher temperature, the system is able to relax in a short 
enough time to a stable single velocity, independent of the 
history of force variation. The behavior does not depend 
significantly on the rate of change of the force since cal- 
culations with different equilibration times give the same 
result as indicated in the figure. The hysteresis observed 
at low temperature is also not affected significantly. How- 
ever, we can not exclude that for large enough equilibra- 
tion times the range of temperature where hysteresis is 
observed is much smaller. The behavior of F a , Ff, and F c 
as a function of temperature is shown in Fig. [?]. Within 
the present model, hysteresis disappears at T ~ Uq < T c 
in contrast with the argumentsB suggesting that smooth 
sliding only occurs above T c . 
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An important feature of the F x Vd curve in Fig. ^ is 
the existence of a "velocity gap ", < v < Vb , which 
leads to stick-and-slip motionB. In the gap region, the 
overlayer cannot sustain a constant sliding velocity and 
stick-and-slip motion is expected when it is driven by a 
spring moving with constant velocity lying in the gap, as 
measured in sliding friction experiments!! Smooth slid- 
ing is possible only for v s > Vb- We have confirmed this 
behavior by calculating the time dependent drift velocity 
and force F(t) when each particle of the overlayer 
is subject to an elastic external force of the form 

F(t) = K s (x cm -v a t) (7) 

where K s is the spring elastic constant and x cm is the 
position of the center of mass of the overlayer. Indeed, 
stick and slip motion is obtained for v s < Vb as shown in 
Fig. U and smooth sliding motion for v s > Vb as shown 
in Fig. g. 

We have also studied the effects of varying the mi- 
croscopic friction parameter 77 on the nonlinear mobility. 
The F xvd characteristics is shown in Fig. for a larger 
value of rj — 0.9 . In this case, the sliding fluid to solid 
transition occur at F c > F a and the return to the pinned 
solid phase at Fb occurs at an effective temperature of the 
overlayer T* below T c , so the results aro, qualitatively 
different from that for rj = 0.6. PerssonH has argued, 
based on hydrodynamical considerations that whenever 
T* < T c is satisfied at Fb , the ratio Fb/F a between the 
kinetic and static friction forces should be a universal 
value ~ 1/2. For the present model, we find instead that 
the ratio Fb/F a depends on 77 and T as shown in Fig. 
[ll]. The ratio approaches 1 for 77 ~ 3 at T = 0.5 . At 
lower temperatures, the results indicate that the ratio 
would asymptotically approach one at larger values of 77. 
This behavior is expected since in the large r\ limit the 
system is overdamped and hysteresis is not expected. In 
fact, we have confirmed through direct calculation that 
when the incrtial term is dropped in the equation of mo- 
tion Eq. ||, hysteresis effects disappear and Fb/F a = 1 
even for the lowest temperature in Fig. [ll]. This should 
be contrasted to the results obtained by PerssonH for the 
Lennard- Jones system where Fb / F a seems to saturate for 
large 77 at w 0.6 for temperatures below T c . 

We should note however that the nature of the high- 
temperature phase in the FK model is quite different 
from the Lennard- Jones system, with power law instead 
of the exponential decay of density correlations as a func- 
tion of separation distance r. This follows from the fixed- 
neighbor constraint in the model of Eq. |l| which ex- 
cludes topological defects, such as dislocations. It is well 
known that thermally excited dislocations can destroy 
crystalline order and lead to exponential decay of density 
correlationsES Therefore, in the absence of dislocations, 
the equilibrium high temperature phase for T > T c and 
F = is a depinned elastic solid and correlations decay as 
a power law. In the sliding state for F > 0, one thus ex- 
pect that when the effective temperature of the overlayer 
T* is larger than T, overheating of the overlayer leads 



to power-law decay of density correlations, whereas long- 
range order should occur when T* = T. This expected 
behavior is confirmed in Fig. |l2] where the finite-size be- 
havior of the structure factor peak for different values of 
F and 77 is shown, corresponding to the different dynami- 
cal phases in Figs. || and [1(| For density correlations that 
decay algebraically as r~ v ° , the normalized peak height 
of S(q) at the commensurate wavevector q — ir should 
scale as a power law S(n)/N 2 cx N~i "° while, it should 
be a constant when long-range order prevailsli-l For the 
fluid phases where T* > T, corresponding to the value 
of F/F = 0.32 in Fig. | and F/F = 0.4 in Fig. [k], 
the loglog plot S(ir) vs N has an apparent linear behav- 
ior with nonzero slope as shown in Fig. |l2|, consistent 
with algebraic order. For the pinned solid phase where 
T* = T at F = 0, the quantity S(n))/N 2 is a constant 
indicating long range order. For the sliding solid phase 
at F/F a = 0.7, the corresponding value for S(tt))/N 2 de- 
creases very slowly with increasing N, indicating that it 
is not a phase with full long range order in contrast with 
the pinned solid phase at F = but much larger system 
sizes would be required to fully confirm this behavior. 

The hydrodynamics arguments used by PerssonH sug- 
gesting a universal value Fb/F a ~ 0.5 are based on the 
existence of a drag force acting on a nucleating domain 
of solid phase surrounded by the fluid phase. Since the 
nature of the fluid phases in the FK and Lennard- Jones 
models are very different, this argument is not valid for 
the FK model and this could account for the different be- 
havior of the ratio Fb/F a . However, in view of the many 
common features shared by the two models as shown in 
Fig. [| and [h], it is also possible that the relation Fb / F a 
~ 0.5 is in fact valid only within a limited range of 77 and 
T, as shown in Fig.[ll], and the Lennard- Jones system 
simply has a broader crossover region than the present 
model. 



VI. CONCLUSIONS 

Comparing our results for the nonlinear response for 
the Frenkel-Kontorova model to .the earlier study based 
on the Lennard-Jones interactiona, we found that all the 
important features such as the appearance of dynamic 
sliding solid and fluid phases, the transition between 
these phases , the strong hysteresis effects, and the ve- 
locity gap for constant velocity sliding are common to 
both models. This indicates that these are robust fea- 
tures that does not depend on the details of the model 
such as many body interactions among the adatoms and 
surface geometry or anisotropies. The strong qualitative 
similarity of the results for mobility in these models to 
the effects observed in macroscopic sliding between two 
surfaces support strongly the idea that the response of 
the lubricant layer to an external force is the key factor 
in controlling boundary lubrication. 

Besides the relevance to the macroscopic sliding fric- 
tion, the dynamic phases generated by the external force 
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are of intrinsic interest and have no counterparts in equi- effects represented by the pinning potential (P. Nozieres 

librium. They are certainly worthy of more detailed and F. Gallet, J. Phys. (Paris) 48, 353 (1987) ). 

study both experimentally and theoretically and maji 
also be relevant to other types of driven lattice systemsEil 
as for vortex lattices in type II superconductors. 
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FIG. 2. Drift velocity Vd and effective temperature T* as a 
function of external force F for r\ — 0.6 and T = 3. Dotted 
lines correspond to free overlayer behavior. 
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FIG. 3. Probability distribution (circles) of particle veloci- 
ties in the stead state for F/F = 0.35 in Fig. m. Solid line is 
a Gaussian fit. 
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FIG. 4. Drift velocity Vd and effective temperature T* as 
a function of external force F for r\ = 0.6 and T = 0.5 and 
periodic boundary conditions. Squares correspond to data for 
increasing F and circles for decreasing F. The arrow indicates 
the minimum velocity Vb- 



FIG. 5. Structure factor S(q) for b/a = 2 and b/a = 1.9 
with r\ — 0.6 and T = 0.5 and free boundary conditions, in 
the sliding (a) and pinned (b) phases. 
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FIG. 6. Drift velocity Vd and effective temperature T* as 
a function of external force F for r\ = 0.6 and T = 1.5. Cir- 
cles and crosses correspond to calculations using equilibration 
times differing by a factor of 10. 



FIG. 8. Time evolution of the Vd and F when the system 
is pulled by a spring at constant velocity v s = 0.01, below Vb, 
for n = 0.6, T = 0.5 and K s = 0.5. 
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FIG. 7. Critical forces F a , Fb and F c as a function of tem- 
perature for rj = 0.6. 
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FIG. 9. Time evolution of the Vd and F when the system 
is pulled by a spring at constant velocity v s = 0.1, above Vb, 
for rj = 0.6, T = 0.5 and K s = 0.5. 



FIG. 10. Drift velocity Vd and effective temperature T* as a 
function of external force F for 17 = 0.9 and T = 0.5. Squares 
correspond to data for increasing F and circles for decreasing 
F. The arrow indicates the minimum velocity v b . 
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FIG. 11. Ratio Fb/F a as a function of r\ and T. 
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FIG. 12. Finite-size behavior of the peak in the structure 
factor S(q) for different values of F in Figs. ^| and |Io| 
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